Preprint typeset in JHBP style - HYPER VERSION 



CMB-Galaxy correlation in Unified Dark Matter 
Scalar Field Cosmologies 



Daniele Bertacca" '*''^, Alvise Raccanelli^, Oliver F. Piattella'^'^ , Davide Pietrobon^, 
Nicola Bartolo"'^, Sabino Matarrese"'^, Tommaso Giannantonio^ '^ 

" Dipartimento di Fisica Galileo Galilei, Universita di Padova , via F. Marzolo, 8 
1-35131 Padova, Italy 

^ INFN Sezione di Padova, via F. Marzolo, 8 1-35131 Padova, Italy 
Institute of Cosmology & Gravitation, University of Portsmouth, Dennis Sciama 

Building, Portsmouth, POl 3FX, United Kingdom 
Departamento de Fisica, Universidade Federal do Espirito Santo, avenida Ferrari 514, 

29075-910 Vitoria, ES, Brazil 

" INFN, sezione di Milano, Via Geloria 16, 20133 Milano, Italy 

^ Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, 
91109 Pasadena CA USA 

^ Excellence Cluster Universe, Technical University Munich, Boltzmannstrafle 2, 
D-85748 Garching bei Miinchen, Germany 

^ Argelander-Institut fur Astronomic dcr Universitdt Bonn, Auf dem Hiigel 71, D-53121 
Bonn, Germany 

E-mails: daniele .bertaccaOpd. inf n. it, alvise . raccanelliOport . ac .uk, 
Oliver . piattellaOgmail . com, davide . pietrobon® jpl . nasa . gov, 
nicola. bartoloOpd. inf n. it, sabino .matarreseOpd. inf n. it, 
tommaso . giaxmantonioOUniverse-cluster . de 

Abstract: We present an analysis of the cross-correlation between the CMB and the 
large-scale structure (LSS) of the Universe in Unified Dark Matter (UDM) scalar field 
cosmologies. We work out the predicted cross-correlation function in UDM models, which 
depends on the speed of sound of the unified component, and compare it with observations 
from six galaxy catalogues (NVSS, HEAO, 2MASS, and SDSS main galaxies, luminous red 
galaxies, and quasars). We sample the value of the speed of sound and perform a likelihood 
analysis, finding that the UDM model is as likely as the ACDM, and is compatible with 
observations for a range of values of Cqo (the value of the sound speed at late times) on 
which structure formation depends. In particular, we obtain an upper bound of < 0.009 
at 95% confidence level, meaning that the ACDM model, for which = 0, is a good fit 
to the data, while the posterior probability distribution peaks at the value = 10~^ . 
Finally, we study the time dependence of the deviation from ACDM via a tomographic 
analysis using a mock redshift distribution and we find that the largest deviation is for 
low-redshift sources, suggesting that future low-z surveys will be best suited to constrain 
UDM models. 
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1. Introduction 

In the last decade, observations of large-scale galaxy distribution [1, 2, 3, 4], the search for 
Type la supernovae (SNIa) [5, 6, 7, 8] and measurements of the Cosmic Microwave Back- 
ground anisotropics (CMB) [9, 10] have been suggesting that two unknown components 
govern the dynamics of the Universe. They are referred to as dark matter (DM), necessary 
to explain the structure formation, and dark energy (DE), that is supposed to drive the 
measured cosmic acceleration (for a review see [11, 12]). However, DM particles have not 
been directly detected yet, although there are hints of their existence [13, 14, 15], and there 
is no theoretically established motivation either for DE [11, 12] or for the tiny cosmological 
constant [16] which would fit cosmological observations [17]. 

A key indication of an accelerated phase in the cosmic history is an excess of power in 
the low multipole region of the CMB angular power spectrum, likely a signatiire of the late 
integrated Sachs- Wolfe (ISW) effect [18]. The late ISW effect is caused by a time-evolving 
gravitational potential and it is interpreted as an effect of a new cosmological component 
which becomes the dominant contribution to the total energy content of the Universe at 
recent epochs. For this reason, it would affect large cosmological scales. Cross-correlating 
the distribution of galaxies with the CMB [19] has been proven to increase the signal-to- 
noise ratio [20, 21] and to be a useful probe of the late-time evolution of the Universe. 

The first correlation between WMAP [22] data and the HEAO-1 X-ray map [23] was 
measured by [20, 21] and then confirmed by the WMAP team [24] using a radio sources 



- 1 - 



map from the NRAO VLA Sky Survey [3]. Further analyses based on different and com- 
plementary techniques strengthened the evidence for the ISW effect [25, 26, 27, 28, 29, 30, 
31, 32, 33, 34, 35] (see [36] for a review) and allowed to test several aspects of cosmological 
models, focusing on properties of DE [37, 38, 39, 40, 41], evolution of radio galaxies and 
their bias [42] and alternatives to general relativity such as DGP [43] and Galileon models 
(sec for example [44, 45, 11]); this technique could in future be applied to even broader 
areas, such as cosmic reionization [46]. 

The standard cosmological model ACDM is a good fit to the current observations, 
but it describes the Universe by means of two unknown components which represent the 
95% of the total energy density. This is theoretically unsatisfactory and alternative models 
which describe DM and DE from an unique prospective have been proposed, such as a single 
component behaving both as dark matter and dark energy, which has been often referred to 
in literature as "Unified Dark Matter" (UDM), or "Quartessence" , see e. g. [47, 48, 49, 50, 
51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72]. Sec [73] for 
a recent review with an up-to-date list of UDM models. In comparison with the standard 
DM -I- DE models (e.g. even the simplest case, with DM and a cosmological constant), 
these models have the advantage that they can describe the dynamics of the Universe with 
a single scalar field which triggers both the accelerated expansion at late times and the 
structure formation at earlier times. Specifically, for these models, we can use Lagrangians 
with a non-canonical kinetic term (see, for example, [74]), namely a Lagrangian which is 
an arbitrary function of the scalar field and of the kinetic term of the scalar field. In [61], 
the authors proposed and studied a technique to reconstruct models where the effective 
speed of sound of the scalar field is small enough that the scalar field can cluster. These 
models avoid the strong time evolution of the gravitational potential and the large ISW 
effect [57] which have been a serious drawback of previously considered models. In this 
reconstruction technique first of all the scalar field Lagrangian was required to be constant 
along the classical trajectories. Specifically, by demanding that L = —A on cosmological 
scales, the background evolution is identical to that of the ACDM [56]. Secondly, from 
this result, they deduced that the energy-momentum tensor of this scalar field is made by 
two components: one behaving like a pressure-less fluid, and the other having negative 
pressure. 

In this paper we consider the cosmological models of Unified Dark Matter introduced in 
[61] focusing on their predictions for the ISW-Large Scale Structure (LSS) cross-correlation. 
We compare them to current observations and constrain the sound speed of the scalar field, 
which is the characterising parameter of the model. 

The paper is organized as follows. In Section 2 we recall the basics of UDM models, 
specifying the Lagrangian and the parametric form of the sound speed for the model we 
consider and its implications for the gravitational potential evolution. In Section 3 we 
describe the theoretical cross-correlation functions for the large scale surveys discussed in 
[33] (NVSS, SDSS main galaxies, LRGs and quasars, HEAO, 2MASS). In Sections 4-5 
we perform the data analysis and discuss the results. Finally in Section 6 we draw our 
conclusions. 
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2. Unified Dark Matter Scalcir field models 



We start recalling the main equations which are useful for the description of most the 
UDM models within the framework of k-essence (see, for example, [74]). Let us consider 
the following action 



S = Sg + S^ = J (fx^ 



(2.1) 



where 

X = -^V^v'V'^^ . (2.2) 

We adopt SttG = = 1 units and the (— , +, +, +) signature for the metric (Greek indices 
run over space-time dimensions, while Latin indices label spatial coordinates). 
The energy-momentum tensor of the scalar field has the following form 

= -^9^ = ^^V,^V„^ + CW,X)g^ . (2.3) 

If X is time-like, describes a perfect fluid and the energy-stress tensor becomes T'^i, = 
{p + p)Ufj,Ui, +pgfj,u which is characterised by the fluid density, p, and pressure, p, once the 
following identifications are made 

p{^,X)=C{v,X) , (2.4) 

p(<^, X) = 2X^^^^ - P{V, X) , (2.5) 

and the four-velocity is defined as 

^- - ^''^ (2.6) 



V2X 

Moreover, let us introduce the equation of state parameter (EoS), defined as the pressure- 
to-density ratio: w = p/ p. 

We will use a flat Priedmann-Lemaitre-Robertson- Walker (FLRW) metric ds^ = —dt'^+ 
a{t)^5ijdx^dx^ = a{r])^{—dr]^ + dijdx^dx^), where a{t) is the scale factor and rj is the confor- 
mal time, and we assume that the energy density of the radiation is negligible at the times 
of interest. Disregarding also the small baryonic component, the background evolution of 
the Universe is completely characterised by the following equations: 

= a'^H^ = laV , (2.7) 

n'-n'^ = a^H = -^a\p + p), (2.8) 

where T-L = a' /a and H = a/a, the dot denoting differentiation w.r.t. the cosmic time t 
whereas a prime w.r.t. the conformal time i]. 

In the background we have that X = 93^/2 = (^'^/(2a^), therefore the equation of 
motion for the homogeneous mode (p{t) becomes 

^ + 2X$^\ (p + ^{3H^) + -^ = 0. (2.9) 

dX dX^J^ dX^ ^' d^dX^ dip ^ ' 
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To the first order in the (scalar) perturbations, we consider small inhomogeneities of the 
scalar field (^(i,x) = ipQ{t) + 6ip{t,x) and we expand the FLRW metric in the longitudinal 
gauge as 

ds'^ = -(1 + 2$)dt2 + (1 _ 2^)a{tf6ijdx'dx^ , (2.10) 

being JT/ = for i ^ j [75]. 

The first-order part of (0 — 0) and (0 — i) component of the energy-stress tensor are 
(see Ref. [76] and Ref. [74]) 



P + P 



+ 3H(p + p)^, (2.11) 



<^r^° = -(p + p)(^) , (2.12) 



where one defines a "speed of sound" Cg as 



dp/dX {dp/dX) + 2X{d^p/dX^ 
Finally, linearizing the Einstein equations one obtains [76, 74] 



2 cin lv%() 

and 



Eqs. (2.14) and (2.15) are sufficient to determine the gravitational potential and the 
perturbation of the scalar field. It is useful to write explicitly the perturbed scalar field as 
a function of the gravitational potential 

^ = 2^^. (2.16) 

Defining two new variables 

where z = a^{p + py/^/{csH), we can recast (2.14) and (2.15) in terms of u and v [74]: 

c,v\=zgy, cv=e(^y (2.18) 

where 6 = l/{csz) = (1 + p/ p)'^/"^ /{^/3a). Starting from (2.18) we arrive at the following 
second order differential equations for u [74]: 

6" 

u" -c^V^u- —u = . (2.19) 
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Unfortunately, wc do not know the exact solution for a generic Lagrangian. However, we 
can consider the asymptotic solutions, i.e. the long-wavelength and the short-wavelength 
perturbations, corresponding to the regimes c^k"^ <C \O"/0\ and c^k"^ ^ l^'V^I) respectively. 
This amounts to consider perturbations on scales much larger or much smaller than the 
effective comoving Jeans length for the gravitational potential [57], which is defined as 
follows: 

X^M = 27r/kj ^ 2t,cI \B/e"\ , (2.20) 

where /cj is the squared Jeans wave number [57]. For a plane wave perturbation u oc 
'Uk(??) exp(zfc • x) in the short-wavelength limit (Cg/c^ S> we obtain from Eq. (2.19) 

- COS (k r cdfi] , (2.21) 

where Ck is a constant of integration at some initial time fj. In the opposite regime, 
neglecting the decaying mode, the long- wavelength solution {c^k^ ^ l^'V^I) 

uy. = Ay,{f})9 f ^, (2.22) 



where is a constant of integration. 



2.1 Matter density and Gravitational potential 

The authors of [61] proposed a technique to construct UDM models where the scalar field 
can have a sound speed small enough to allow for structure formation and to avoid a strong 
integrated Sachs- Wolfe effect [18] in the CMB anisotropics which typically plagues UDM 
models [50, 51, 57] (sec also [64, 69]). By applying this technique we can explicitly derive 
the functional forms of the energy density and gravitational potential for the class of model 
we are interested in. 

2.1.1 Energy density 

In particular, starting from the following scalar field Lagrangian £, [61] 



X) = cosh (if) /l- — 

' (l-c2,)sinh(7<^)[l + (l-c2^)sinh2(7(^)]V A 

A [(l-c^)2sinh2(7<^) + l-cy ^2 2^^ 



(l-cl) l + (l-cysinh2(<^) 

1 /9 

where 7 = [3/[4(l - c^)]] ' , and once the initial value of (p is fixed at early times, 
the scalar field Lagrangian is constant and equal to —A along the classical trajectories. 
Specifically, when C = p = —A on cosmological scales, the background is identical to the 
one of the ACDM. Indeed, if we consider the equation of motion of the scalar field, and if 
p = —A, we easily obtain 

p [a{t)] = pDM(a = l)a"^ + A = poM + PA , (2.24) 
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where p\ behaves hke a cosmological constant (pA = const.) and poM behaves hke a DM 
component (pDM oc a~^). This result imphes that we can think of the stress-tensor of the 
scalar field as being made of two components: one behaving like a pressureless fluid, and 
the other having negative pressure. Therefore the integration constant /9DM(a = 1) can be 
interpreted as the "dark matter" component today with Jlmo = Pdm{o, = l)/(3iJo^) and 
^Ao = pa/ (3i?o^) the density parameters of DM and DE. In this case, the parametric form 
for the sound speed is [61] 

CsHa) = ^ , ,f , (2.25) 



^^AO + (1 - Coo^)^^mOa" 



where Coo is the value of the sound speed when a — > oo. Let us emphasise that when a ^ 0, 
then Cg — 7- 0. 

Specifically, our class of UDM models allows the value w = —1 for a oo. In 
other words, they admit an effective cosmological constant energy density at late times. 
Therefore, in order to compare the predictions of our UDM model with observational data, 
we follow similarly the prescription used in [67] , where the density contrast of the clustering 
fluid is (5dm = Sp/pBM and where pdm = P — PA, by definition, is the only component 
of the scalar field density which clusters (for adiabatic EoS see [62, 67]). According to 
this notation we can derive the expression for the density contrast as function of the 
gravitational potential for scales smaller than the cosmological horizon and z < z^ec, where 
Zrec is the recombination redshift (zrec ~ 10^): 

A ^P(.^;z) _ fc^^(k;z) 

<^DM (k, .) - ^^^^^^^ ^ - ^^^^^ ^ . (2.26) 

2.1.2 Gravitational potential 

Let us stress here an important point which holds for the models that we are considering 
(i.e. those which reproduce the same expansion history of the ACDM Universe): the 
gravitational potential evolves in the same way as in a ACDM Universe for those modes 
with wavelengths larger than the sound horizon. In fact, we can write the usual solution 
[77, 78] 

$fc«i/,,(k; z) = (l - I' a\fj)dn^ , (2.27) 

for k"^ <^ 1/Aj(r/); here i?k — (9/10)<I>p(k)Tjn(fc), where $p(k) is the primordial gravita- 
tional potential at large scales, set during inflation, sec e.g. [78], and T^{k) is the matter 
transfer function suggested in [79] , that includes baryons (a similar approach is been used 
also in Ref. [69]). 

In this regime, k"^ < 1/Xj{r]), we can write $(k;z) ~ {9/10)^p(}s.)T^(k){l + z)D{z), 
where D{z) = g{z)/{l + z) is the growth factor and g{z) the growth suppression factor for 
a ACDM Universe. A very good approximation for g{z) as a function of redshift z is given 
in Refs. [80, 81, 82]: 



1+^ 



-^_^^^a(^) 



70 



(2.28) 
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withj^,, = n,,,o{i+zf/E^{z), = nAo/EHz), E{z) = {i+z)n{z)/no = [n^oii + + ^^ao]'^' 

and r^mOi ^Ao = 1 ~ ^mOi the prcscnt-day density parameters of non-relativistic matter 
and cosmological constant, respectively. We have normalised the growth suppression factor 
so that g = I early on when — > 1 and Q,\ — > 0. 
In the opposite regime we have 

1/2 



^fe>i/Aj(k;r/) 



1 



p + p 



(j])Ck{'n)cos (^k J_ Cs{f})dr]^ , when A;^ » 1/Aj(?7) . 

(2.29) 



In Eq. (2.29) we have Ck(^) {9/10)^p{k)T„^{k)C , where C* is a constant of integration 



1 



C = 2^ " ,^ . (2.30) 

The value of C is obtained under the approximation that for rj < fj one can use the long 
wavelength solution (2.27). In other words, we are assuming that there is an epoch in the 
past when the sound speed is very close to zero (a condition that is necessary in unified 
models to allow for structure formation). Notice that Eq. (2.29) clearly shows that the 
gravitational potential is oscillating and decaying in time for the perturbations inside the 
sound horizon. Moreover, through Eqs. (2.26), we can draw the same conclusion also for 
5dm (k,z). 

In order to connect these two regimes, we can introduce a scale-dependent growth factor 
D{k;z) as $(k;z) ~ (9/10)$p(k)rm(fe)(l + z)D{k;z). In the long-wavelength regime, 
k^ ^ 1/Aj(77), D{k;z) just reduces to D{z) = g{z)/{l + z), according to the previous 
discussion. In this way one can naturally define a transfer function, Tudm, as 

Dik; z) = Tudm(A;; ri{z))D{z) , (2.31) 

or, in other words, 

$(k; r/) = ruDM(A;; ?7)$fe«i/Aj (k; r/) . (2.32) 

In particular, we impose that TuDM(fe)??) = 1 for jy < r^rec, where r/rec is some epoch when 
the Universe is matter dominated and the radiation is negligible (usually later than the 
recombination epoch). Before this epoch the sound speed must be very close to zero to 
allow for structure formation. 

In Fig. 1 we show the evolution of Tudm as a function of the scale factor for = 0.1 h 
Mpc~^ and of the wavenumber for a = 1. In order to investigate their behaviour at 
different scales, we study the evolution of the growth factor and the gravitational potential 
for some values of A;. In Fig. 2 we compare the ACDM case with the UDM one with 

= 10^^, 10~^, 10~^; for the smallest value of c^, the behaviour coincides with the one 
of the ACDM model, and the difference from ACDM increases with k (for the gravitational 
potential plots sec also Rcf. [64]). These plots show that for large values of Cqo (when Jeans 
length becomes very large), the structure formation is affected and becomes unstable, and 
this will cause a decrease in the cross-correlation. 



-7- 



/.5 p • • r 
1.0 

'\ 

\ 

0.5 - 
0.0 - 



ACDM 

cj=10' 

cj=10' 

cJ=10' 



\ 



— 0.5 I ^ , , \ , ^ ^ \ ^ ^ ^ \ ^ ^ ^ \ ^ ^ ^ 

0.0 0.2 0.4 0.6 0.8 1.0 

a 



y g I ' ' ' I ' ' ' I ' ' ' I 

ACDM 

cj=10-^ 

cj=10-'' 

cj=10"^ 



1.0 



II \ 
a 

5;- 0.5- \ 
d 



0.0 - \ ,/ 



/ 



-0.5 



0.00 0.02 0.04 0.06 0.08 0.10 

k [h Mpcr^] 



Figure 1: Illustrative plot of Tudm as a function of the scale factor for k = 0.1 h Mpc^^ (top 
panel) and as a function of the wavenumber for a — 1 (bottom panel). For each panel, from left to 
right, the curves correspond to the choice ~ 10^^, 10^'^, 10^^ and ACDM model. 



2.1.3 Analytic approximation 



We found a simple analytical expression for the UDM transfer function which fits very well 
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Figure 2: Growth factor D-k{z) (left column) and normalized gravitational potential ^v{z) ~ 
<i>(k; z)/$(0; 10^) (right column) as function of redshift for different Fourier modes. Top to bottom, 
the sound speed has been set to = 10"^, 10"^, 10"''. For k < 10"^ h Mpc"^ the UDM follows 
the ACDM behaviour (black solid line). For larger modes the discrepancy increases with the sound 
speed. 



the numerical results given by (see also [72]) 

TjjY,M{k-r^) = jo{A{ri)k) , (2.33) 

where 

n da 

Aiv)= cM)d,' = c^ Ha\/l + il-c- ).a-3 ' ^^'^'^ 

where v = Omo/^Ao ~ 0.387 and = HqQ\o{1 + ua'^). Since ^(r/) is not analytical 

and we are going to consider < 10"'^ [57], it is convenient to expand Cg in Eq. (2.25) in 



Taylor series near Cqq = 0, i.e. 



; , + \ , '^'^'^',3/2 ^^ + O {cl) , (2.35) 



where u = i^mo/^AO ~ 0.387. The approximation 

Coo 



(2.36) 



is therefore very good, being the truncation error of order c^. Plugging Eq. (2.36) in 
Eq. (2.34), changing the integration variable to the scale factor a and choosing a = as 
lower integration limit we find 



Aia) Coo £ „ , r. = I da^:^ . (2.37) 



da _ Coo r ^ 



The integration can be performed analytically and the result is: 

Hon^^Aja) VStt 1 a2 _ + ^2/3 ^ ( y/^2a-v^ 

^ "'18^ + 6^^'' (a + ^V3)2 +^^^^^^^(^-3- ^1/3 J • 

(2.38) 

In order to check the accuracy of the approximation of our fitting function (2.33), let 
us define the following relation: 

f 2 \ _ l-^Fit ~ -^Numeric I _ |[/?UDM(A;;a)dA; da] - [/($/$fc«i/A,)(A;; a)dA; da]\ 



I °° -^Numeric [/($/$fe^i/Aj ) (A); a)dA; da] 

(2.39) 

Prom Pig. 3 we can notice that at most we have a relative error about 3% for = lO^^ 
Now, using the transfer function of Eq. (2.33), let us analyse in detail the behaviour 

of the gravitational potential, the growth factor and the power spectrum of (5dm- 

In Pig. 4 we show the power spectrum of the UDM energy density component that 

clusters: 



Pik;z) = -T^{k)T^^uik;z)l-^] I — ) P^{k), (2.40) 



where 

'mO \ , -3 



\D{z = 0)7 \Hi 
is the primordial 3D power spectrum of the potential field, 

($p(k)$p(k')) = (27r)3<5(k + k')P$(A:) , (2.42) 

6h is the expression given by Eq. (A3) of Ref. [82] and A = 1.06 is computed normalizing 
the power spectrum to ag = 0.78 (see [31]). 
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Figure 3: Percentage difference between numerical solution and the analytic approximation used 
in our analysis {AI/I, Eq. (2.39)) as function of c^. It can be seen that the analytic approximation 
is very good and the discrepancy with respect to the exact solution is at most few percent. 



Fig. 4 confirms that, when k > kj, $(k;z), D{k;z) and P{k;z) oscillate and decay 
both as a function of time and k. As a result, the general trend is that the possible 
appearance of a sound speed significantly different from zero at late times for these UDM 
models corresponds to the appearance of a Jeans' length under which the dark fluid does 
not cluster anymore. At small scales, only if the sound speed is small enough, UDM 
reproduces ACDM (see also [57, 64, 69]). 

3. Cross correlation between the CMB and the Large Scale Structure 

The integrated Sachs- Wolfe effect is an energy shift which a photon exhibits when it crosses 
a gravitational potential evolving in time. Since in a Universe made of pressure-less mat- 
ter the gravitational potential is constant, in the standard cosmological scenario we may 
distinguish two effects: one at early time, early ISW effect, during the transition from a 
radiation-dominated Universe to a matter dominated one, and one at late time, late ISW 
effect when the evolution of the Universe starts feeling the presence of a cosmological con- 
stant. The late ISW takes place at recent epoch and affects large cosmological scales: this 
reflects into a bump at low multipoles in the angular power spectrum of the CMB. It is 
well known that this region has a large intrinsic uncertainty due to the cosmic variance and 
the constraining power results pretty low. The cosmological constant affects the structure 
formation at large scales as well: cross-correlating the distribution of galaxies with the 
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Figure 4: Power spectrum P(k) of the UDM models for a wide range of 



CMB [19] has been proven to increase the signal-to-noise ratio [20, 21] and to be a useful 
probe of the late time evolution of the Universe. 

In particular this is true in the case of UDM models. Indeed it is important to stress 
that in the UDM models there are two simple but important aspects: first, the fiuid which 
triggers the accelerated expansion at late times is also the one which has to cluster in order 
to produce the structures we see today. Second, from the last scattering to the present 
epoch, the energy density of the Universe is dominated by a single dark fluid, and therefore 
the gravitational potential evolution is determined by the background and perturbation 
evolution of just such a fluid. As a result, the possible appearance of a sound speed 
significantly different from zero at late times corresponds to the appearance of a Jeans 
length (or a sound horizon) under which the dark fluid does not cluster anymore, causing 
a strong evolution in time of the gravitational potential and of the fractional overdensity 
Sdm which start to oscillate and decay. Thus, besides having the possibility of a strong 
ISW effect for UDM models [57], it is interesting to analyse the possible changes of the 
cross-correlation between the distribution of galaxies and CMB signal in UDM models due 
to a non-negligible sound speed, for example, with respect to the ACDM model. 

Following Refs. [37] and [38], consider a held x{'y) as a function of the angular position 
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7 on the sky to be a weighted projection of the potential field $(x; z): 

Sx{j) = J dzW^f,z)Hf,z), (3.1) 

where the weight can also, in general, include differential operators acting on the field. 
The angular cross-correlation between two observed fields x and x' is given by 

c-'(^) = {6x{^)dx'{n = E^^(^ ■ ^o^c-r' • (3-2) 

The brackets () denote averaging over all direction 7 and 7', satisfying the condition 7-7' = 
cos(i?), Pe is the Legendre polynomial and Cf^' is the angular cross power spectrum, which 
is given by 

where the weight If{k) is given by 

Itik) = ^T^{k) j dz{l + z)TuDM(fc; z)D{z)W''{k- z)u{kx{z)) , (3.4) 

and je{kxiz)) are the spherical Bessel functions. Here x(z) = dz/H{a) is the 
comoving distance to redshift z. 

When cross-correlating CMB and LSS to measure the ISW effect, Sx = ST = @, where 
is the CMB temperature fluctuation, and 5x' = (5„g(7, z), where Sng = Sug/ug is the 
galaxy overdensity observed and where ng{^, z) is the number of sources in the sky which 
traces the matter distribution [24]. The window functions read 

W^ik- z) = -2rcMB e-(^) ^lM(l + ^)TuDM(fc;.)i^(.)] 

az 

(3/2) ^golo(l + .) ^- 

where e"'^^^) is the visibility function, which accounts for the effect of reionization; h{z) is 
the galaxy bias and dN{z)/dz is the mean number of galaxies per steradian with redshift 
z within dz. Therefore, the integrand functions lf{k) and I^''{k) are respectively: 

rT/,N 9 [ -r(^) d [(1 + ^;)ruDM(fc; z)D{z)\ . 

h {k) = -^TcMB Tn^[k) dz e ^^^^ — Jl[kx{z)] , (3.7) 

Ii'ik) = J dz b{z)^^T^^M{k;z)D{z)ji[kx{z)\ . (3.8) 



Finally, the Legendre transform of the cross power spectrum " is the angular cross- 
correlation function, which can also be defined as the average 



".(,?) = (e(7i)<5n,(72)) = fo(e(7l)5„^(72)), (3.9) 

where the CMB temperature fluctuation and the galaxy overdensity 6ng are observed in 
any directions 71, 72 separated by an angle i9. For simplicity, we have assumed that the 
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galaxy bias b is constant for each data set. The theoretical cross-correlation functions for 
the large scale surveys discussed in [33] (NVSS, SDSS main galaxies, LRGs, and quasars, 
HEAO, 2MASS) are plotted in Fig. 5. The actual shape of each catalogue depends on 
the redshift distribution of the sources, as shown in top panels both for the two-point 
correlation function and its Legendre transformation, however in all cases considered, the 
UDM model with < 10"^ is extremely close to the ACDM one. In the right lower 
panel, the particular case of SDSS main galaxies catalogue is computed for different values 
of the parameter = 10~^, . . . , 10~^ to highlight the effect of the speed of sound on the 
CCF. The bottom left panel displays the redshift distributions of the catalogues to help 
understanding their properties. 



ACDM 



ACDM 




S 10 15 20 

B [deg] 

Survey Redshift Distribution 




2 MASS 

SDSS gal 

LRG 

NVSS 

HEAO 

QSO 



0.30 

o.zs 

0.20 
0. 15 
0.10 
0.05 

0.00 







3MASS 
SDSS pai 
LRG 
■ NVSS 



HEAO 

QSO 



100 

I 



Theoretical CCF for SDSS 




5 10 IS 

g [deg] 



Figure 5: Top: cross-correlation function between WMAP temperature CMB map and 6 cata- 
logues of sources both in real space (left panel) and harmonic space (right panel). The redshift 
distribution of each catalogue is shown in the bottom-left panel whereas the effect of is shown 
for the specific case of the SDSS galaxy survey (bottom-right). 



The ISW effect is sensitive to the time derivative of the gravitational potential, hence 
for the UDM model we discuss, it is expected to be higher for large values of the sound 
speed: the cross-correlation signal in these cases is lower than in the ACDM model. This 
effect can be explained through the oscillatory behavior of the density contrast and of the 
gravitational potential when k > kj (see Eqs. (2.20), (2.21) and (2.26)). Indeed this is due 
to two possible factors: i) through the amplitude of the density contrast that decays over 
time: the presence of a non-negligible speed of sound prevents the structure formation on 
scale smaller than the Jeans' length (see Eq. 2.20), so that one expects a cross-correlation 



- 14 - 



in these cases lower than in the ACDM model; ii) the rapid fluctuations between positive 
and negative values produce a temporal mismatch between $' and 5dm which determines 
a lowering of the cross-correlation signal. 

In the next section, we perform a analysis to determine which value of the sound 
speed fits better the data. The measured 2-point correlation functions as well as the 
covariance matrix between the catalogues are publicly available^, whereas the theoretical 
curves are computed implementing the algorithm described in [31] following the specifics 
in [33]. This analysis is not meant to fully span the parameter space of the model, which is 
beyond the purpose of this paper. We assume the WMAP best-fit flat cosmology [83] for 
Ob^^, ^ch^, Ho, Us, As and J7a and vary Cqo only, which is a peculiar feature of the UDM 
model under investigation. 



4. Data analysis 

We use the ISW data from [33] , which were obtained cross-correlating six galaxy catalogues 
with the WMAP maps of the CMB. This data set probes the distribution of the large-scale 
structure in a wide redshift range < z < 2, using different parts of the electro- magnetic 
spectrum. The used catalogues arc: the 2MASS infrared survey, the optical SDSS, which 
includes main galaxies, luminous red galaxies and quasars, the radio-galaxies from NVSS 
and the X-ray background of HE AO. Although not all the redshift distributions arc known 
with a great accuracy, this represents a first approximation towards a true tomographic 
reconstruction of the potentials. 

All maps were coarsely pixellated on the sphere with a resolution of 0.9 deg. Then the 
real-space cross-correlation functions (CCFs) were measured in angular bins of 1 deg for 
Odeg < < 12 deg, so that the full data set for the six catalogues consists of 78 points 

The full, strongly non-diagonal, covariance matrix Cij representing the uncertainties 
on these points was calcolated in three ways, using a model-independent jack-knife (JK) 
approach and two different Monte Carlo methods. In the first (MCI), the real LSS maps 
were cross-correlated with 5000 random CMB maps generated using Gaussian random 
seeds. In the second {MC2), also the LSS maps were randomly produced, based on their 
redshift distributions. These methods give comparable results, although it was found that 
the JK approach is rather unreliable, due to strong dependencies on the particular followed 
procedure. Here we decided to use the MCI, in order to conservatively avoid having to use 
the redshift distributions of the catalogues in the calculation. 

We then calculate the likelihood of a given model by calculating its theoretical CCF 

[cj^] , based on the measured redshift distributions of each catalogue and on linear 
\ /theo ^ 

theory. The likelihood is then obtained as £ oc where 

- E [(cf-),, - (/"-U [C-]., [(?'•)„,. - (cf-) J . (4.1) 



^http://www.usm.um-muenchen.de/~ tommaso/iswdata.php 



We assume that the galaxy bias is constant for each catalogue and we allow it to vary 
around the value 6fid reported by previous authors and confirmed by [33], in the range 
0.5 6fid < b < 2 6fid, and marginalising over this additional degree of freedom. We discuss 
our findings in the next section. 

5. Results and discussion 

The total likelihood for the speed of sound is shown in Fig. 6. The 95% confidence level 
limits are € [0., 0.009]. This reveals that despite the large number of surveys combined 
together the constraining power is weak and we hit the priors. The ACDM model is a 
good fit to the data, but the distribution peaks to the value = 10~^. We check the 
stability of the solution against the bias correction repeating the analysis assuming the 
fiducial value for each catalogue. While the different treatment of the uncertainty on how 
baryon trace cold dark matter marginally affects the confidence region at 68%, the 95% 
bounds are unchanged. 

400 I J 




0.000 0.002 0.004 0.006 0.008 0.010 

2 
C ^ 

Figure 6: Total likelihood of the parameter c^, in the UDM niodel resulting from the combined 
analysis of the publicly available galaxy surveys. The horizontal dot-dashed line marks the 95% 
confidence interval. 

In order to investigate whether this peak is a signature of any peculiar process hap- 
pening at a specific redshift we performed the analysis for each catalogue separately. The 
result is plotted in Fig. 7. For many surveys the likelihood results monotonically decreasing 
with the speed of sound, peaking at the value which recovers the standard model. Two 
of them show actually a peak at = 10"'^: HEAO and LRG. By looking at the redshift 
distribution of the sources. Fig. 5, no common feature emerges as a clear indication of 
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some physical process. The two distributions peak roughly at z ~ 0.6, but the LRG one is 
extremely sharp, whereas HE AO very broad. 
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Figure 7: Constraints on set by each catalogue. Only two catalogues, HEAO and LRG, show 
a mildly peaked distribution at ~ 0.001. The solid line shows the likelihood marginalised over 
the nuisance bias parameter, whereas the dashed line is obtained fixing the bias of each catalogue 
to its nominal value. 



We conclude then that current galaxy surveys are not very sensitive to the speed of 
sound of the class of models we studied, the main reason being a not accurate enough 
redshift characterisation of the samples, which would allow us to perform a tomographic 
study of the time evolution of the cross-correlation function. 

5.1 Tomographic analysis 

In order to test for which range in redshift we could have the best opportunity to discrim- 
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inate between UDM models and ACDM, we perform a tomographic analysis, using five 
mock redshift distributions, following the model in Ref. [37] (see the left panel of Fig. 8). 
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Figure 8: Top-left: Mock redshift distribution and redshift bins used to analyse the redshift 
dependence of ISW. Top-right, bottom-left and bottom-right: angular (cross-)power spectrum for 
each redshift bin for the ACDM model, the choice = fO""* and the relative difference. As 
expected, the most sensitive bin is the first one, peaked at z ~ 0.4; interestingly the largest difference 
between the models occurs at £ ^ f 00. 



Since the standard model coincides with the choice = 0, with = 10~^ being 
already a good approximation, and a UDM model with = 10~^ seems to be slightly 
favoured by current observations, we compare the two cross-correlation functions computed 
for the same redshift slices as a function of the peak of the galaxy distribution, Zp, looking 
for a signature of the different time evolution of the two models. In the right panel of 
Fig. 8 are our results; as one can see, and as we expected (see Ref. [64]), the maximum of 
the sensibility is for low-redshift surveys, so future tests of this model should be performed 
with tracers at low redshift. Indeed, sources at lower Zp emit light that strongly feels the 
decay and the oscillations of the Newtonian potential, because this potential is sensitive to 
(-2^5 Coo); that decreases with time, and to the presence of an effective f^A, that plays the 
role of DE [64, 57]. 
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6. Conclusions 

In this work we have considered the cosmological models of UDM introduced in [61] focusing 
on predictions for the two-point correlation function between CMB and galaxy distribution. 
The main goal of this work is to investigate the ISW effect within the framework of UDM 
models based on a scalar field and understand which constraints we can put on the sound 
speed of the scalar field, crucial parameter of these models, using the CMB-LSS correlation. 
Analysing the cross-correlation between different source catalogs (NVSS, SDSS, HEAO, 
2MASS) and the WMAP CMB temperature map as discussed in [33], we studied its changes 
due to a non-negligible sound speed and compared them to the standard ACDM scenario. 

As expected, wc found that for a sound speed smaller than 10^'^, the UDM model gives 
the same prediction of the ACDM one, whilst for > 10"*^ the model starts to behave in 
a different way. In particular, this effect can be explained through the oscillatory behaviour 
of the density contrast and of the gravitational potential when k > kj (see Eqs. (2.20), 
(2.21) and (2.26)). Indeed this is due to two factors: i) through the amplitude of the 
density contrast that decays over time: the presence of a non negligible speed of sound 
prevents the structure formation on scale smaller than the Jeans' length (see Eq. 2.20), 
so that one expects a cross-correlation in these cases lower than in the ACDM model; ii) 
the rapid fluctuations between positive and negative values produce a temporal mismatch 
between and 5dm which determines a lowering of the cross-correlation signal. 

In conclusion, the comparison of theoretical predictions with real data has shown that 
current surveys are not precise enough to decisively distinguish between specific classes 
of late-time evolution of the Universe. Future surveys will have a better redshift charac- 
terisation of sources so that they will allow a tomographic study. We investigated this 
prospective creating a mock catalogue with a realistic redshift distribution. In particular, 
from tomographic analysis, the maximum of the sensibility is for low-redshift surveys. In- 
deed, sources at lower Zp emit light that strongly feels the decay and the oscillations of the 
Newtonian potential, because this potential is sensitive to kj{z,Coo), that decreases with 
time, and to the presence of an effective Oa, that plays the role of DE. 
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